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Abstract 

We give the general analytic solutions derived from the low energy string effective action for 
four dimensional Friedmann-Robertson- Walker models with dilaton and antisymmetric tensor field, 
considering both long and short wavelength modes of the i?-field. The presence of a homogeneous 
H-tie\d significantly modifies the evolution of the scale factor and dilaton. In particular it places a 
lower bound on the allowed value of the dilaton. The scale factor also has a lower bound but our 
solutions remain singular as they all contain regions where the spacetime curvature diverges signalling 
a breakdown in the validity of the effective action. We extend our results to the simplest Bianchi I 
metric in higher dimensions with only two scale factors. We again give the general analytic solutions 
for long and short wavelength modes for the H field restricted to the three dimensional space, which 
produces an anisotropic expansion. In the case of H field radiation (wavelengths within the Hubble 
length) we obtain the usual four dimensional radiation dominated FRW model as the unique late 
time attractor. 



"edmundjc@central.susx.ac.uk 
''a.lahiri@central.susx.ac.uk 
'^dwandsQstar. susx.ac.uk 



1 Introduction 



String inspired cosmology is currently attracting a great deal of attention. The most favored starting 
point in any analysis is the low energy string effective action from which the lowest order string beta- 
function equations can be derived Q. These equations, for the closed string, consist of three long range 
fields, the dilaton (j), the Kalb-Ramond field strength Hfj_y\, and the graviton, all arising out of the 
massless excitation of the string. In addition there is a constant related to the central charge of the string 
theory which vanishes in the critical number of dimensions 10 or 26. The fact that only the massless 
excited state is used suggests that the effective action is not a valid description for probing the highest 
energies associated with string theory. However, we may hope that through the beta-function equations 
we are investigating physics associated with events from say the string scale down to the GUT scale. 
Such an approach has already been adopted by a number of authors [||, ^, ^, ^, ^, |^, |^ 0, [T2| . 

Veneziano and his collaborators have emphasized the possible importance of the duality symmetry 
which characterizes the equations of string cosmology [||, 11, |l^. They point out the possibility that 



inflation can occur without relying on a potential energy density. Rather, the duality transforms, which 
relate the dilaton and the metric, lead to a possible inflationary mechanism. In [^] it is claimed that 
the inclusion of the ff-field does not seem to change the underlying properties of these duality-related 
cosmologies, although they also point out possible problems with exiting inflation in these models. 
Cosmological solutions with a dilaton and a non-trivial if-field have been obtained by TseytlinB, for 
curved maximally symmetric spaces in string theory, with a non-zero central charge deficit. (See j^] for 
details and complete references). 



In [12 1, the authors study the outcome of string-dominated cosmology in a four dimensional Friedmann- 
Robertson- Walker (FRW) spacetime including a homogeneous i?-field {Hfj_y\) as well as a non-zero crit- 
ical charge deficit V (some of the particular solutions were previously obtained by Tseytlin |^]). Their 
results are based on a phase-plane analysis. In this paper, we present the general analytic solutions, 
including both long and short wavelength solutions for the //-field, but setting V = and then show 
how this can be extended to simple anisotropic models in higher dimensions. 

In section 2 we introduce the low energy equations of motion in the string frame and demonstrate how 
they can be related via simple conformal transforms to equations in other frames. Section 3 concentrates 
on four dimensional FRW models and describes the complete solution for a homogeneous //-field, 
reproducing where applicable previous solutions found in the literature |P, |l2|. We go on to describe 
radiation solutions in which the //-field has a spatial dependence on small scales (i.e. wavelengths much 
smaller than the Hubble length). Because the metric (and the dilaton) is homogeneous and isotropic 
we require the //-field energy-momentum tensor to be homogeneous and isotropic on average, and 
demonstrate how at late times we recover the general relativistic results for radiation and curvature 
dominated models. In section 4 we consider the influence of homogeneous and radiation solutions of 
the //-field va. D = A + n cosmologies where the metric tensor is decomposed into the direct product 
of a four dimensional FRW metric and an n dimensional metric. We give general analytic solutions for 
the simplest version of the D dimensional Bianchi I models where there are just two scale factors. The 
results indicate how the //-field can produce an anisotropic expansion leading to one scale expanding 
whilst the other contracts. Finally in section 5 we summarize our main results. 

2 String action 

We shall take as our starting point the low energy /^-dimensional string effective action [Kt] 



S=^ I d^xV^e-"^ R + {V(l)f -V - -^//2 



2«;|, 



12 



+ / d^X^/^£matter , (2.1) 
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where k^-, = SttGd and we adopt the sign conventions denoted (+++) by Misner, Thorne and Wheeler 
1 14]. <p is the dilaton field determining the strength of the gravitational coupling and = H^i^xH^^^^ 
where H^^x = d]^^B^x] ■ 

The variation of this action with respect to the g^y^ B^j^y and (j), respectively, yields the field equations 



j^;: {y^r + [gX'' - ax'') ^av.^ , (2.2) 

(e-'^Rf"'^) = , (2.3) 



2n(P + R-{V(j))^ -V -^H"^ = 0, (2.4) 

where T^i, is the energy-momentum tensor derived from the matter Lagrangian. 

The effect of certain types of "stringy matter" has been considered elsewhere in the literature |11|. 
Specific schemes of compactification, not to mention the chosen gauge symmetries of the theory, will 
determine the behavior (and number) of both bosonic and fermionic matter fields. It is not inconceivable 
that in some 'matter-dominated era' of the stringy epoch of the universe these matter fields will play 
a part in determining the cosmological evolution. In particular, the antisymmetric tensor may be 
considered a matter field, and as we shall see it plays a significant role in string cosmology. In favor 
of considering the possible effects of this field on the cosmology we shall ignore all other contributions 
from the matter Lagrangian. 

The charge deficit ^ is a constant proportional to D — 26 for the bosonic string and D — 10 in the 
heterotic or superstring. We will set V = in our analysis. This may well be necessary for a consistent 
theory either by choosing the appropriate number of spatial dimensions or due to cancellation with 
contributions from other matter fields. Even for non-zero V this should be an increasingly good approx- 
imation at early times if the curvature and/or kinetic energy densities are large, V <^ R, {V (j))'^ , H'^ , but 
we would need to consider its effect at late times in an expanding universe. 

2.1 Conformal frames 

These field equations are similar to those found in Brans-Dicke gravity []l5| with the Brans-Dicke pa- 
rameter io = —1. This is only strictly true in the absence of the H-&eld as in Brans-Dicke gravity it 
is assumed that the energy-momentum tensor of all fields (other than the Brans-Dicke field, <I> = e~'^) 
are minimally coupled to the metric g^y. While the energy-momentum tensor of other matter fields 
are assumed to be conserved with respect to this metric (the string metric), so that 'Vf^T^'^ = 0, we 
cannot define an energy-momentum tensor solely in terms of the i7-field and the string metric which 
is conserved independently of the dilaton. This is just a consequence of the equation of motion for H 
(Eq. (^I^)) which has an explicit dependence upon e"^. 

Hfj.ux is only minimally coupled in the conformally related metric 

g^i, = exp ( ^'^ ) g^i, , (2.5) 



which we might call the B-metric, in which we findEI 

y^^M^A ^ Q ^2.6) 



^Note that the 3- form //^^a has a conformally invariant definition in terms of the potential 
form has indices raised by a particular metric. 



2 



Another particularly useful metric to introduce is the Einstein metric[16| 

20 



g^u = exp 



D-2 



(2.7) 



In this frame the action appears simply as the Einstein-Hilbert action of general relativity in D- 
dimensions, while the dilaton appears simply as a matter field, albeit one interacting with the other 
matter fields, 
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The corresponding field equations are then those for interacting fields in general relativity; 
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where the terms on the right-hand side of the Einstein equations correspond to 

D(j) 



rpV 



4 



4 (^)r; 



^2 {V)fu 
I^D -'-fj, 



exp 
1 



D-2 



1 



D-2 

1 

— exp 
12 ^ 



1, 



v~\k 



-,9^9 



--Fexp 



40 
D-2 
20 
D-2 



2 



9, 



(2.9) 
(2.10) 
(2.11) 
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(2.15) 



the energy-momentum tensors for the matter, dilaton and iJ-fields and potential V respectively. The 
total energy-momentum must be conserved of course by the Ricci identity, but there are interactions 
between these four components. Henceforth, as remarked earlier, we shall set T^u and V to be zero. 



3 Isotropic D=4 solutions 

Firstly we will consider the behavior of 4-dimensional homogeneous and isotropic cosmologies for which 
the most general metric is the Friedmann-Robertson- Walker metric 

ds^ = -dt^ + a\t)(^^ + r^dnA , (3.1) 
\ 1 — Kr^ / 

= a^(r?)(-dr,^ + ^+rW) , (3.2) 

in terms of the conformally invariant time coordinate, r], with k = +1,0, —1 for spatially closed, flat or 
open models respectively. Just as we take the metric to be homogeneous we shall also assume that the 
dilaton has no spatial dependence, = 0(ry). 
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With D = 4 the conformal transform to the Einstein frame ffives a rescaled scale factor 
The metric field equations are simplest in terms of the Einstein metric where we have the familiar 
constraint equation (the 00-component of Eq. ( |2.9D with ^ = and = 0) 

- 3 f 5 + 4) = ^^^To" + ^-^^^o") , (3.3) 

where prime denotes differentiation with respect to r]. 

In 4-dimensions the i/-field equation of motion, Eq. (|2.10|), is solved by the Ansatz, 



where is the antisymmetric 4-form (obeying Vp?'"'^^'' = 0) and the integrability condition, dy^H^xi^-^ - 
0, becomes the new equation of motion for h 

U h + 2V''(l)V^h = . (3.5) 

Thus h evolves as a massless scalar field coupled to the dilaton (except in the B-frame where □ h = 0). 
The same interaction appears in the dilaton equation of motion (Eq. ( |2.11| ), with V = 0) as 



n 6 = e^'t> 



{Vhf . (3.6) 
Thus we have two interacting scalar fields whose energy-momentum tensors are given by 

^'^'^^t; = l{glr''-\r,t'')<t^,x<t^,., (3.7) 

^^'^"^T; = ]^[g>;;g'^--]^r,t^y^^hxh,. (3.8) 

These equations simplify considerably in two cases. 
3.1 Homogeneous solution h{ri) 

Thus far we have allowed for the possibility of a spatial dependence of the //-field. However as we 
have already restricted ourselves to considering a homogeneous metric and dilaton field this will only 
be consistent with choosing a source that is at least homogeneous on average. Indeed as far as we are 
aware, the only case that has been considered to date|12] is that of a strictly isotropic ff-field, where 
h = h{r]), and thus Hq^^^ = 0. 

In this case the equation of motion for the //-field (Eq. (p.5|)) becomes 



a 



h" + ( 2^ + 2</.' ) /i' = , (3.9) 



which is easily integrated to give e^'^TPh' = itL where L is a non-negative constant. It is this kinetic 
energy density of the //-field that drives the dilaton. Note that back in the string frame this solution 
corresponds to = ^1? jt^ . Thus it will dominate any charge deficit V in the dilaton equation of 
motion, Eq. ( |2.4| ), as a ^ 0. 

Because both h and (/> are functions only of time, we can define a new scalar field '0(^) where 
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The energy-momentum tensor in the Einstein frame is then simply that for this single minimally coupled 
field 

The equation of motion for a homogeneous minimally coupled field is 

V'" + 2^V' = 0, (3.12) 
a 

which can be integrated to give a^V' = where X is a positive constant. (We will not consider the 
trivial case where K = {).) Thus in the Einstein frame we have an isotropic perfect stiff fluid whose 
energy density K^p = -k^{ WTgO + ^""if^) = K'^/Aa^. 



The constraint Eq. (|3.3| ) for the Einstein scale factor, a, can then be integrated to give [117 1 



where we define 



I tan(r/ — ?7o)| for ^ = +1 , 
'T'i'n) = { for A; = , (3-14) 

tanh(?7 — r/o)| for k = —1 . 

We emphasize that in the Einstein frame the scale factor evolves in a wholly unremarkable fashion ||lO|| . 
We have a singularity at ry = r/o with 2 = and the models expand away from it for r] > rjQ or collapse 
towards it for rj < rjQ. Only closed models can turn around, and these recollapse at = r^o ±vr/2. There 
are no bounce solutions. Notice also that the behavior of the Einstein scale factor is independent of L, 
and thus is the same in vacuum, i.e. without the presence of the -ff-field, as it is with the i?-field. 
Combining the first integrals for -0 and the //-field with the definition of dip we also have 



- - . (3.15) 



a4 



This too can be integrated to give 



±v/3 

for L = , 



2K 



where r=K is an integration constant. Note that the evolution in the presence of the antisymmetric tensor 
field is quite distinct from the vacuum behavior. In particular there is a lower bound on the dilaton, 
g2(/> > K^JL?. By contrast, there are two distinct branches in vacuum where the dilaton is either 
monotonically increasing or decreasing. 

We can use this to recover the scale factor in the string frame 



Kp. [i + kT"^)-^ forL = 0, 

' Li^ +(^) 1 + A:r2 ^ for L / . 



,2 _ J 

Lti, 

2^3 

The evolution of the scale factor and dilaton in different cases is shown in Figs. (1-4). The singular 
behavior of the conformal factor, e"^ , at the initial singularity in the Einstein frame produces cosmologies 
which have a minimum non-zero value for the scale factor in the string frame. In the presence of the 
//-field, a diverges both as r/ — > ±oo (or as r/ — > ryo i 7r/2 for A; > 0) and as — > r/o, so all models 
"bounce" , although they are still singular in the sense that the Ricci curvature diverges. 
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t -> 



t -> 



Fig. 1: The decelerated branches (left) and the accelerated branches of the scale factor a (solid line) and 
dilaton e"^ (dotted) in a spatially flat FRW universe with vanishing (or constant) h. 

Note that the vacuum solutions again exhibit two distinct branches corresponding to (j) and a either 
monotonically increasing or decreasing. The k = models correspond simply to power-law solutions 
shown in Fig. 1. We can write the solutions in terms of the proper time in the string frame by integrating 
dt = adrj. The two solutions are then 

• e'f' (x\t- tol^^""^ and a (x \t - to\+^/^; 

• e<^ oc |t - to\-^-^ and a (x \t - to\'^'^/^. 

These two branches, corresponding to a decelerated or accelerated scale factor, are related by the duality 
transformation 1 18 1 a — > 1/a and — > e"^ 



fa 
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For L 7^ the A: = solution corresponds to the power-law vacuum solutions at early and late 
times and we find that it smoothly interpolates between the accelerated and decelerated branches. (See 
Fig. 2.) In some sense then it could be described as a "self-dual" solution. 

• as ?? — > 7/0) we have e*^ oc |t — to\^^^^ and a oc |t — tol"""^^^ as t — > to; 

• as 1] ^ ibcxD, we have e"^ oc Itl^^"*" and a oc Itl-*-/^ as t — > itcxo. 

This is precisely the behavior found by Goldwirth and Perry [p!^ in their phase-plane analysis. The 
vacuum solutions correspond to particular solutions (here corresponding to the limiting behavior where 
r^, is either infinite or zero) found previously Q. 

The solutions in spatially curved models approach the flat space results only near r/Q. At late times 
r ^ 1 in open (k = —1) models (Fig. 3) and thus the dilaton becomes frozen-in at a fixed value as the 
curvature dominates the evolution and and we approach the Einstein result. In closed models (Fig. 4) 
where r — > oo as r/ — r/o — > ±7r/2, the scale factor diverges in a finite proper time. Thus although these 
models bounce, and therefore must undergo a period of inflation (a > 0) in the string frame, they still 
become curvature dominated at late times. 



6 



t -> 



Fig. 2: Scale factor a (solid line) and dilaton e'* (dotted) in a spatially flat FRW universe with homogeneous h. 




t -> 



Fig. 3: Scale factor a (solid line) and dilaton e"^ (dotted) in an open (fc = —1) FRW universe with homogeneous 
h. 
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t -> 



Fig. 4: Scale factor a (solid line) and dilaton e*^ (dotted) in a closed (fc = +1) FRW universe with homogeneous 
h. 

3.2 Radiation solution, (V/i)^ — 

It is also possible to consider cases where the //-field does have a spatial dependence. Because our 
metric (and dilaton) is homogeneous and isotropic we will require that the i/-field energy-momentum 
tensor is also homogeneous and isotropic on average. It is natural (in fiat space) to decompose any field 
into Fourier modes, h = J2qhq{r]) ^wi^QiX^) where Qi is a spatial comoving three-vector. We see then 
that the preceding case corresponds to the long wavelength mode where = Qi ~^ 0- The other 
case in which we can solve the equation of motion is where g — > oo where effects of spacetime curvature 
can be neglected and the usual Minkowski spacetime result holds. 

Specifically we find that, in the B-frame, for S> d" /a, k/a^, the equation of motion for hq{ri) 
reduces to 

(dhq)" + q\dhq) . (3.18) 

Thus, for a single short wavelength mode, in the limit q oo, we have hg = rriqe^'^'^ /d where rUq is a 
constant. This corresponds to an energy-momentum tensor for the i/-field in the Einstein frame 

(^)t; = q.q^"^ , (3.19) 

where q^ is a null 4-vector with q^ = q. Clearly this is not isotropic for a single wave-vector, g^, but if 
we consider an isotropic distribution of wave- vectors we find 

-K^{ (^^f^) = Sn'^i ^^'^fi) (no sum), (3.20) 

- f. (-) 

the usual result for radiation in a FRW universe, where M = J rn^q^dq. 

Note that the fluid is trace-free, and thus conformally invariant. This means that the energy- 
momentum of the radiation is conserved in all conformal frames and so the i/-field is decoupled from 
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the dilaton which appears as a minimany coupled scalar field in the Einstein frame. Thus the first 
integral of its equation of motion, for a homogeneous (pirj), gives 



A 



(3.22) 



Just like the ip field for homogeneous h{r]), the energy-momentum of the <f) field behaves like a stiff fluid, 
with 



— K 



2 WfO 



^2 (^)fi 

ASP ' 



(no sum) 



(3.23) 
(3.24) 



Once again the constraint Eq. 



be integrated|17] to give 



now for two non-interacting fluids, one radiation, one stiff, can 
_n K + Mt t 



V3 1 + A:r2 



(3.25) 



using the time coordinate deflned in Eq. (3.14). (As in the homogeneous case, this is just the familiar 
behavior for a FRW universe in general relativity with matter obeying the strong energy condition 
and thus singular for all models at r 



0.) This in turn allows the equation for (j) ■, Eq. ( |3.22| ), to be 
integrated, yielding and thus, via the conformal transformation, the scale factor in the string frame. 



±V3 



K + Mt t 



^/3 1 + A;r2 ' 



where we have introduced 



s{ri) 



MT{ri) 



(3.26) 
(3.27) 

(3.28) 



K + MT{r]) 

Notice how the evolution of the dilaton fleld is now similar to that in the vacuum case except it is 
determined by the function 5(77) rather than T{rj). At early times (?7 ~ %) s is proportional to r but, 
unlike the vacuum case, s ^ 1 as r ^ 00 and so the fleld becomes frozen in at late times in the flat 
model as well as the open model (where r ^ 1 and thus s — > M/{K + M)). Thus we recover the late 
time general relativistic results for radiation and curvature dominated models respectively. 

In common with the vacuum solutions, there are two distinct branches (Fig. 5) with the dilaton 
monotonically increasing from zero (the decelerated branch) or decreasing from inflnity (the accelerated 
branch) at rj = r]Q. For Mr <C |i^| the evolution is essentially identical to that in the vacuum case with 
two branches where a — > cxd as r ^ when (/> — > 00, but a — > when — > 0. Thus we have no "bounce" 
solution interpolating between the vacuum branches. 
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t -> 



Fig. 5: The decelerated branches (left) and the accelerated branches (right) of the scale factor a (solid line) 
and dilaton e"^ (dotted) in a spatially flat FRW universe with short wavelength h. 



4 Anisotropic D = n + 4 solutions 

Having investigated the homogeneous and radiation solutions for the H-&eld in four dimensions, we 
return to Eqs. (|2.8|)~(|2.14|) in the Einstein frame where we will consider a metric tensor in D = 4 + n 



-( 2.14 ) in the Einstein frame where we will consider a metric tensor in D 
dimensions that can be decomposed into the direct product form 



dp + a'^{t)dxidx'' + h\t)dxjdx-' 



(4.1) 



where we let i run from 1 to 3 and J run from 4 to n + 3. The scale factors a and h thus refer to the 
3-space and n-space split respectively. We have chosen these spaces to be spatially flat, thus this is a 
Bianchi I metric. The procedure we outline here could also be applied to the general D-dimensional 
Bianchi I metric, but we restrict our discussion here to a two-scale factor model to avoid introducing 



too many degrees of freedom. Eq. (4.1) is of course related to the original string metric through the 
conformal transformation Eq. ( |2.7| ). 

The Einstein evolution equations, Eq. (p.S|), for the two scale factors written in terms of d = 
{dd/dt)/a and $ = {db/dt)/b are then 



d + {3d + n/?)d 
/3+ (3d + n/3)/j 
plus the constraint equation 



n + 2 
n + 2 



(p+ (n - l)p-nq) , 
{p + 2q-3p) , 



9 ■ n(n — 1) A9 9 ~ 

3d^ + 3nd$ + — = njjp 



(4.2) 
(4.3) 

(4.4) 



where p = —Tq, p = T- and q = fj (no sum over i or J). Note that any isotropic stiff fluid for which 
p = p = q makes no contribution to the right-hand-side of the evolution equations, entering only into 
the constraint equation. 
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We will adopt the simplest extension of our = 4 Ansatz for the iJ-field, 



where 

e'^^^'^ = -^^5^^51515^ . (4.6) 

Note that whereas this Ansatz included all the solutions in 4 dimensions, this represents only one of 
many degrees of freedom for the ff-field in the D = 4 + n case. An antisymmetric field in the higher- 
dimensional space corresponds to a whole range of fields in the 4-dimensional reduced theory|19|. At a 



classical level we can obtain self-consistent results if these other fields are set to zero initially, though in 
a full quantum treatment we would have to consider precisely how the massive modes of these fields are 
excited in the extra dimensions. Our choice of Ansatz makes any dependence of the function h on the 
n-space coordinates, x"^, irrelevant as this cannot affect H^^'^^ and so we need consider only h = 
A similar approach of demanding the //-field live only in three space dimensions was considered in Q 
who found numerically that in the Einstein frame the effect was to drive that space to a large size while 
the other spaces remained of order the Planck scale. 

Using (|4.5|) the integrability condition becomes the equation of motion for h{t,x^); 



3 

li,u=0 



f^^df, (s36-"e-^<^/("+2);j^^^ ^ . (4.7) 



We can decompose the field h into its Fourier components, h{t, x*) = hq{t) exp(iqjX*) so that the equation 
of motion for the homogeneous function hq is 



hq+ I 3d -nf3+ —-4)]hq + -^hq = . (4.8) 



The final equation of motion is that for the dilaton in the Einstein frame, Eq.( |2.1l| ), driven by 
^2 ^ _6exp(8(/)/(n + 2))h-^'"' {V hf giving 

. . g4</./(n+2) 2 

(/) + (3d + n(5)4) = ^ — (V/ij . (4.9) 

Note that having chosen all the fields to be independent of the coordinates on the n-space we could 
replace the D dimensional action by an effective theory written in terms of the four dimensional part 
of the metric in the string frame 



where the effective D = A dilaton is given by 



i?4 + (V99)'-n(V(ln6))2-li/2 



(4.10) 



= ^ ' (4-11) 



and the scale factors of the extra dimensions just act as massless (moduli) fields. However to emphasize 
the dynamical evolution of the n-space we shall treat this scale factor on an equal footing with that of 
the 3-space. Note that the effective Z) = 4 Einstein frame would not be the same as the D = 4 + n 
dimensional Einstein frame due to the modified dilaton. 
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4.1 Homogeneous solution h(t) 

The equation of motion for h for long-wavelength modes (where — > 0) becomes 

h+(3a-nP + ^-(j)^h = 0, (4.12) 



or 



n + 2 



where L is a positive constant. From Eq. ( 2.14 ) and making use of Eqs. ( [4. 51) and ( 4.13| ) we have 

H^XkH""^^ = if^or i/ = {0,4,5,..n + 3} , 



21^ . 

HanH^'"' = -^H fori audi = {1,2, 3} 



Hence in Eq. ( p. 14 ) we obtain 



^^"'^fl ^ 4ph = ^e-^<^/('^+2) (no sum), (4.15) 
Kl^^)fj ^ 'cW = -^e-^<^/("+^) (no sum), (4.16) 



which means that the -ff-field acts like an anisotropic fluid satisfying 

PH =PH = -QH ■ (4.17) 

Note that, just as in the D = A case, in the string frame is inversely proportional to the square of 
the volume of the 3-space {H^ oc a~^). 

There is also a contribution to the energy- momentum from the dilaton field given by Eq. ( 2.13| ), 

. = (4.18) 

Kl^^lf = kJ, Wfy = ^Jl^ (no sum), (4.19) 

which means that the (p field acts like an isotropic stiff fluid, 

p<t> = P<l> = q<p , (4.20) 

and thus, as remarked in the preceding section, drops out of the evolution equations for a and 

We can substitute the energy- momentum tensors into Eqs. ( |2.11| ), (4.2) and ( f4.3D to obtain the 
following equations of motion, 

<^ + (3a + n0)<}) = ^e-^'^/("+2) , 
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while the constraint equation becomes 



+ 



2(n + 2) 4S6 



-A(t>/(n+2) 



(4.21) 



Notice how the presence of the //-field on the right-hand-side of these evolution equations tends to 
drive (f) and a in a positive direction, but drives j3 negative. Thus it produces shear and an anisotropic 
expansion. 

Introducing a new time coordinate ^ through 



we obtain first integrals for the equations of motion, 



a 







2(n + 2) 

I j2^i<j>/{n+2) 



(n + 2) 



62" 



(4.22) 

(4.23) 
(4.24) 
(4.25) 



and 



h = ±L{i + ih). (4.26) 

where cj)' = dcjy/d^ etc, and ^^-.iai^p and are constants of integration. In fact at least one of these 
is redundant as the origin of the variable ^ is clearly arbitrary as ^ is only defined by the differential 
relation in Eq.( ^.22| ). Henceforth we shall take = so that /?' = at ^ = 0. We can solve for /3' by 
differentiating Eq. ( [4.25|) and substituting in for (j)' from Eq. (|4.23| ) to obtain 







(4.27) 



h (n + 2)e + + C ' 

where C is another integration constant related to the others via the constraint Eq. ( 4.21| ). 

(n + 2)C7 = 2(n + 2)^^ - -in^il . (4.28) 

There is another important constraint which emerges. From Eqs. ( [4.25| ) and ( 4.27 ), by demanding 
that the combination 52ng-40/{n+2) j-gj^^^ins non-negative we obtain 



n 



+ 2)e + + C < , 



(4.29) 



which means that the allowed range of ^ is bounded by m_ < ^ < m+, where m± are the roots of the 
above expression; 

2^0 ^ A 



m± 



(n + 2) n + 2 ' 

and we have introduced A = [Sn^^^ _ 2^^2ji/2 Q^gg^j-^y solutions only exist for 3n^2 > 2^,. 



(4.30) 



We solve Eq. ([4.27|) to obtain 



where is a constant and 



KO = ^o(e-"^~)^-(m+-e) 



9+ 



5± 



2(n + 2) 



1 + 



(4.31) 
(4.32) 
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We obtain similar solutions for a and 4> using Eqs. ( [4 .241 ), ( |4.23D and ( 4.31| ) 

a{i) = ao(e-m_f-(m+-er , 
= e^o -m^y (m+ 

where clq and (/)o are constants (with 6^e~^'^°/^""'"'^^ = L ), p± and s± being given by 



n 



P± 



4(n + 2) 



It 



A 



1 ± 



A 



(4.33) 
(4.34) 



(4.35) 
(4.36) 



It is clear from Eqs. ( 4.31| ), (|4.33 ) and ( 4.34| ) that the behavior of a, b and e'^ is similar in each case, 
the differences emerging in the exponents of the (^ — m±). We need to know how these solutions appear 
in the string frame as this is where the original theory has emerged from. This is straightforward to do. 
Using Eqs. (|2.7|) we see that the scale factors in the string frame are given by 



a(0 

where the constants u± and v± are 



ao(e-m_)«-(m+-e)"+ , 



(4.37) 
(4.38) 



u± 



v± 



1± 



A 



2A 



Thus the qualitative behavior of the solutions is quite simple, approaching power laws as ^ — > m± with 
the exponents u± and v±. 

Using Eqs. ( |2.7| ) and ( [4.22 ) we see that the time (dt) in the string frame is related to (d^) by 

a^h-^'e'f'di , 



dt 



where 



oc (C 



w± 



-l+W- 



[rriA 



-1+10+ 



di, 



(4.39) 
(4.40) 



1 + 3ti-|- — nv± + s± 



1 ± 



3n^„ 



4A 



(4.41) 



Although we cannot in general integrate this relation to obtain t{^) in closed form, we can solve for the 
time in the string frame in the limits ^ m± to give 



t oc 1^ — m± 



(4.42) 



We see that ^ — > m± in finite proper time if w± > 0, which means that ^ — > m_|_ in finite proper time if 
< 0, and rri- in a proper time if > 0. Thus the proper time interval is always semi-infinite. 

(The case ^q, = is excluded by the requirement that A 7^ 0.) 

FVom Eqs. (|1|), (IISTI) , (g^H) and (|I3|) we have the limiting behavior of the solutions as power 

laws 



a{t) oc 
b{t) oc 

P0W 



oc 



(4.43) 
(4.44) 
(4.45) 
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Thus we need to know the behavior of the ratio u±/w± in order to understand the late time behavior 
of the solutions for a{t). Assuming that we are dealing with events at late proper time as ^ ^ rrij^ and 
t — > oo (i.e. > 0) then we have 

^ = 4^ , (4.46) 

which is bounded by 

1 u+ 1 

77 < — < ^ • 4.47 
3 ^3 

The upper bound corresponds to the case C,/) = where h =constant and we recover the late-time 
behavior of the isotropic D = 4 solutions. This implies that in the string frame we do not obtain 
inflationary solutions (by which we mean cPa/dt'^ > 0) at late times. This of course does not imply that 
we do not obtain inflationary solutions at earlier times. On the contrary, all the solutions which initially 
contract {da/dt < 0) but then expand at late times {da/dt > 0) must undergo a period of accelerated 
expansion. 



Comparison with the isotropic D = 4 solutions discussed in Section (3^) is easiest if we use the 



conformally invariant time coordinate defined by drj = dt/a, which is related to ^ via 

(4.48) 

where ?/* — ??o = Oo/n^aL. Note that for < the coordinate r] runs — oo to tjq and for > we must 
have rjo < r] < oo, i.e. the range for tj is always semi- infinite, coinciding with the semi-infinite range for 



7] 


-m 




rj* 


-m 





the proper time. To avoid any ambiguity we introduce the non- negative variable r defined in Eq. ( 3.14 ) 
but here restricted to the spatially flat {k = 0) case 

r = |r/-r/o|, (4.49) 
which runs from oo to for ^q. < and from to cxd for > 0. We then have 



T 




(4.50) 

(4.51) 
(4.52) 

Note that 6 is a monotonic function of r. However the qualitative behavior of a depends on the 
value of n^a/^<f> and can be separated into two cases. 

1. A < In^ol (requiring \nS^a\ ^ \S,<f>\)- '^i'^) is monotonically increasing function; 

2. A > \n^a\ (requiring \n^^a\ > \C<f>\)- contracts initially, "bounces" and then expands. 

The latter case includes the D = 4 isotropic case (shown in Fig. 2) where = and b remains constant. 
We then have A/n^a = \/3 and we recover the solutions given in Eqs. (l3l^ ) & (IsTtI) independent of 



n, the number of extra dimensions while they remain static. Another example is when we start with 
an isotropic contraction (a/a = b/b < 0) at r = as shown in Fig. 6. The size of the n-space continues 
to decrease whereas the -ff-field prevents the 3-dimensional space collapsing and it then expands at late 
times. 



15 



t -> 



Fig. 6: Scale factors a (solid), h (dashed) and dilaton e"^ (dotted) in a = 4 + 6 Bianchi I universe with 
homogeneous h, starting from an approximately isotropic state a.i t — Q. 

The critical case where A = \nE,a\ = corresponds to the case where the scale factor a remains 
finite but non-zero at r = 0. Thus the curvature of the 4-dimensional spacetime is nonsingular, although 
it must be emphasized that both the dilaton and n-space scale factor remain singular. This "wormhole" 
solution in the D = b case has been discussed recently by Behrndt and Forste[^]. 

Although the (4 + n)-dimensional dilaton (p displays the same range of qualitative behavior as the 



scale factor a, the behavior of the effective 4-dimensional dilation ip defined in Eg. ( 4. 11 ) is much more 
restrictive, evolving as 

( / T \ / T \ ' 



Thus it always has a minimum value at r = r*. As in the isotropic D = A case, the presence of the ho- 
mogeneous antisymmetric tensor field introduces a minimum allowed value for the effective gravitational 
coupling constant. 

4.2 Radiation solutions, (V/i)^ = 

As in the D = A isotropic case we can seek solutions corresponding to the short wavelength limit of the 
Fourier decomposition of h[ri,x^) = hq{r]) exp(zgjX*), where it is again convenient to write the equation 
of motion for a Fourier mode hq using the conformally invariant time coordinate, defined such that 
drj = dt/a = dt/a. We then have in the extreme short wavelength limit, g — > oo, 

and thus _ 

where niq is a constant and is a null 4- vector restricted to the 4-dimensional metric d's^ = o?{—d'rf + 
dx^). 
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The corresponding energy-momentum tensor in the Einstein frame is 

which when averaged over an isotropic distribution (with respect to the 3-space) gives the anisotropic 
(with respect to the whole (n + 3)-space) radiation fluid 

^ -l~PH = :^^, (4.57) 

4(Wf;^) ^ 4p^ = i-^ (no sum), (4.58) 

«D ^"^T-j = K%qH = (no sum) . (4.59) 
The dilaton field again acts as an isotropic stiff fluid with 

which does not appear in the Einstein evolution equations. 

Substituting this ff-field solution into the evolution Eqs. ([4.9|), ([4.2|) & (|4.3|) we have 



^ + (3d + n/j) = , (4.61) 
d + (3d + n/3) d = , (4.62) 



3a46" 

/3 + (3d + n/j) /3 = , (4.63) 



while the constraint Eq. (4.4) becomes 



3d2 + 3nd/3 + = + . (4.64) 

2 2(n + 2) ^ ^ 

We can again obtain first integrals of all three evolution equations this time using the conformally 
invariant time, r/, so that 

,/ 1 ^ . N 

= i ('? + '?»). ("C) 

/3' = (4.67) 

where 77,^, ija and ry^ are constants of integration. 

Using the same technique as in the homogeneous case, differentiating Eq.( [4.66 ) and substituting in 
for /?' using Eq. (|]6|) we obtain a second order equation for a{r]) whose first integral gives 

r/2 + (2r]a + ^nr]j3)ri + C 
The constant C is given in terms of the other constants from the constraint Eq. ( [1.64D as 

C = rg+ Snri^r,^ + '-7]p - ^^^^^vl ■ (4.69) 
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Comparing Eq. ( [4.68| ) with Eq. ( [4.66| ) we see that 

~2^„ ^ ^ (^^2 ^ (2r/„ + 3nr]p)r] + c) , (4.70) 

For each choice of integration constants we have two distinct solutions. We have two semi-infinite 
intervals; one where t] approaches m_ from —oo, and the other for rj > m_|_ approaching +cxd. m± are 
the roots of the above expression; 

3 

'm± = -r]a - -nr]p ± A , (4.71) 

and 



A = i^3n(n + 2)r?2 + -l-r?2 (472) 

Integrating the first-order equation for a(r/), Eq.( 4.6^ ), and similar expressions for j3 and finally 
yields 

a{rj) = aQ\r] — m-\^~\rj — mj^\^+ ^ (4.73) 
6(77) = bQ\ri — m-\'^~\rj — mj^\'^+ , (4.74) 
= e'^o|ry-m_|^-|??-m+|*+ , (4.75) 



with the exponents 



= (4,78) 

These are related to the solutions for the two scale factors in the original string frame by the 
conformal transformation, Eq. (^]^), which gives 

a{ri) = aol?? - m_|"-|?7 - , (4.79) 

h{ri) = bQ\r]-m^\'"-\r]-m+\'"+ , (4.80) 

(4.81) 



with the exponents 



- .5^. (4.83, 



2A 



(4.84) 



As in the homogeneous case the conformal time rj cannot in general be given in closed form as a 
function of the proper time in the string frame, and can be given instead only by the differential relation 
dt = adrj. However, unlike the homogeneous case, we can identify a unique late/early time behavior as 
ry — > ±00. We then find dt oc rj^-'^'^+drj and thus |t| oc r/^ giving 

a oc \r]\ (X , (4.85) 

b, e'^ constant . (4.86) 
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Thus the unique late time attractor for rj > rrij^ is the famihar solution for an expanding D = A radiation 
dominated universe. For t] < rri- the asymptotic solution at early times represents the usual solution 
for a collapsing radiation dominated universe. Both these solutions are singular as r/ — > m± with a — > 
for M-t > or a ^ oo for u± < 0. 

To compare these solutions with the isotropic case discussed in Section ( |3.2| ) we can write them in 
terms of the non- negative time variable introduced for the spatially flat D = A metric in Eq. ( 3.14| ); 



^ ^ I m_ - 7? for r]<m^ , 
I ry — for 7] > m-|_ . 

When rj < rri-, r decreases as the solutions approach the singularity at r = 0, while when rj > m_|_, r 



increases away from the singularity. Then s{r]), defined in Eq. (3.2g|), is given by 



and we can write 



siv) - ^ , (4.88) 



a = aorV2(^ + A)V25(2A-nB)/2^ ^^gg) 
b = bos^+^ , (4.90) 



where 



B ^ ^ . (4.93) 

SO that (n + 2)(2A^ + nB'^) = 3. We can now identify the isotropic case with b remaining constant and 
A = —B = ±\/^/{n + 2), which yields the results of Section (3.2). 



5 Conclusions 

We have written down the general solutions to the low-energy equations of motion derived from string 
theory for a four dimensional Friedmann-Robertson- Walker cosmology including a homogeneous dilaton 
field and antisymmetric tensor field. 

In the absence of the anti-symmetric tensor field {H^yx = 0) we recover the known vacuum results 
for the evolution of the dilaton. In spatially flat models we have two distinct branches; the decelerated 
branch where both the scale factor, a, and the dilaton (j) grow monotonically from zero at the curvature 
singularity, and the accelerated branch where a and (p decrease from infinity at the singularity. These 
solutions are invoked in the so-called "pre-big-bang" cosmologies based on the accelerated branch ex- 
panding from a and (p zero with low curvature, and approaching the high curvature regime, where it 
must evolve into the expanding decelerated branch with a and (p ^ oo as the curvature again vanishes. 
But achieving such a smooth transition in the context of the low energy eff'ective theory has proved 
difficult even with the inclusion of a dilaton potential|13]. 

We have shown that the presence of the homogeneous i?-field introduces a minimum value for both 
the FRW scale factor and the dilaton in L) = 4 models. Spatially flat models have a ^ oo both at 
the spacetime curvature singularity at ?? = r/o and in the low curvature limit (?? ±00). We no longer 
find two distinct branches but rather a solution which smoothly interpolates between the accelerated 
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vacuum branch at high curvatures and the decelerated branch at low curvatures. Our analytic solutions 
confirm some of the results of the phase-plane analysis by Goldwirth and Perry Thus the extreme 
weak coupling limit, (/> — > 0, is attained only when H^^x = 0. This is in contrast to the results of |11| 



because our i^-field does not obey the requirement for 0{d,d) symmetry in their solutions in that our 
B^i, is not homogeneous. Rather, we consider solutions where it is the field H^^x = d^^Bj^x] which is 
homogeneous. These solutions are in fact related to the homogeneous vacuum solutions by an SL(2, R) 
duality transformation [^] where the dilaton and the H-G.eld are considered as the real and imaginary 
parts of a single complex fielda. 

We find qualitatively similar solutions when we consider 4-|-n dimensions. We give analytic solutions 
for a Bianchi I model with two scale factors where the antisymmetric tensor field acts as a scalar field 
on only three of the spatial dimensions. It tends to accelerate their expansion with respect to the other 
n dimensions producing an anisotropic expansion even from isotropic initial conditions. The effective 
four dimensional dilaton always has a minimum value. We recover the isotropic Z) = 4 solutions in the 
limit that the second scale factor b =constant. 

By treating the antisymmetric tensor field as a scalar field on the four dimensional sub-space we can 
consider the short wavelength modes which act like radiation as well as the long wavelength homogeneous 
modes. We have also given analytic solutions in the case where this radiation is isotropic on three spatial 
dimensions. The dynamical effect of the ff-field is then negligible near the curvature singularity and 
we see both the accelerated and decelerated branches seen in vacuum. However we find a unique late 
time attractor solution in spatially flat models corresponding to the usual D = 4 radiation dominated 
solution in general relativity with a tx t^/^ and — ^constant (and b — ^constant in the two scale factor 
anisotropic model). The dividing line between the short and long wavelength regimes corresponds to 
wavelengths within or outside the Hubble length. Although this is a frame dependent quantity, for 
power-law evolution, a tx |?7 — ?7o|^) (which we find in all asymptotic limits) any given wavelength must 
become "long" as r/ — > 779 and conversely any wavelength eventually becomes "short" as |r/| — > 00. 

In models with non-zero spatial curvature (considered only in the isotropic D = 4 case) we find 
that although our solutions approach the fiat solutions near the spacetime curvature singularity they 
become dominated by the spatial curvature as r/ ^ ±00 where the spacetime curvature becomes small. 
This is hardly surprising if one considers the evolution of the metric in the Einstein frame, of which we 
have made extensive use throughout, where these string models correspond simply to a universe with a 
massless scalar field (plus radiation for the short wavelength ff-field) and so there can be no inflation 
in this frame. 

It is important to emphasize the limited cosmological era in which the results presented here may 
be valid. While the appearance of a minimum value for the scale factor of three dimensional space is 
intriguing, the string metric is still in general singular at rj = r]Q. By solving the fleld equations only 
to lowest order in the string coupling constant a' we are neglecting terms of order a' R^^xkR^^^^^ with 
respect to terms such as i?, H^, (V(/))^, etc. in the fleld equations. As all our solutions approach power- 
law evolution (cj oc |?/ — ^I^') as ry — > 770, then these terms inevitably become divergent, so our solutions 
will only be good approximations to the true evolution when the spacetime curvature is sufflciently 
small. On the other hand we have neglected any potential for the dilaton, assuming V <C (Vi;^)^. If we 
hope to recover the standard hot big bang cosmology at late times then we will have to include matter 
flelds which presumably provide a potential to flx the present day value of the dilaton (and thus satisfy 
observational limits on the allowed variation of the gravitational coupling strength within the solar 
system today). Indeed if this potential for the dilaton or for any other fields produces an infiationary 
era any memory of a preceding stringy era would be all but erased. Neglecting the dilaton potential is 
only likely to be valid while kinetic terms dominate. As we find oc a^^ (for homogeneous modes) 
this will only to be valid at sufficiently small a. 



We are particularly grateful to N. Kaloper for drawing this to our attention. 
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In summary, we have given general analytic solutions for the evolution of an early, but sufficiently 
low-energy, stringy era where the massless bosonic fields dominate the dynamics. We have shown that 
the presence of the antisymmetric tensor field has a dramatic effect on the evolution of the dilaton in four 
dimensions and can also produce an anisotropic expansion in higher dimensional models preferentially 
expanding three spatial dimensions. 
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